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Abstract 

We examine the effects of possible phenomenological interactions between dark energy and dark matter on cosmological 
parameters and their efficiency in solving the coincidence problem. We work with two simple parameterizations of the dynamical 
dark energy equation of state and the constant dark energy equation of state. Using observational data coming from the new 182 
Gold type la supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson 
Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, 
we perform a statistical joint analysis of different forms of phenomenological interaction between dark energy and dark matter. 
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I. INTRODUCTION 



Our universe is believed undergoing an accelerated expansion driven by a yet unknown dark energy (DE) [1]- [•")]. 
Much effort has been devoted to understand the nature and the origin of the DE. The leading interpretation of such 
a DE is a cosmological constant with equation of state (EoS) w = —1. Although this interpretation is consistent 
with observational data, at the fundamental level it fails to be convincing. The vacuum energy density falls far below 
the value predicted by any sensible quantum field theory, and it unavoidably leads to the coincidence problem, i.e., 
"why are the vacuum and matter energy densities of precisely the same order today?" . To overcome the coincidence 
problem, some sophisticated dynamical DE models either relating the DE to a scalar field called quintessence with 
w > —1, or to an exotic field called phantom with w < ~1 [(>] have been put forward to replace the cosmological 
constant. But it is doubtful that there is a clear winner in sight to explain the DE and solve the coincidence problem. 

DE contributes a significant fraction of the content of the universe, it is natural to consider its interaction with the 
remaining fields of the Standard Model in the framework of field theory. The possibility that DE and DM can interact 
has been studied in [7]~[1G], among others. It has been shown that the coupling between a DE (or quintessence) field 
and DM can provide a mechanism to alleviate the coincidence problem [7, 10]. A suitable choice of the coupling, 
motivated by holographic arguments, can also lead to the crossing of the phantom barrier which separates models 
with EoS w > —1 from models with w < —I [11, 12]. In addition, it has been argued that an appropriate interaction 
between DE and DM can influence the perturbation dynamics and affect the lowest multipoles of the CMB spectrum 
[13, 17]. Arguments using structure formation of galaxies suggested that the strength of the coupling could be as 
large as the QED fine structure constant [13, 18]. More recently, it was shown that such an interaction could be 
inferred from the expansion history of the universe, as manifested in, e.g., the supernova data together with CMB and 
large-scale structure[l!)]. Nevertheless, the observational limits on the strength of such an interaction remain weak 
[20] . In addition, it was suggested that the dynamical equilibrium of collapsed structures would be affected by the 
coupling of DE to DM [21, 22]. The basic idea is that the virial theorem is distorted by the non-conservation of mass 
caused by the coupling. This problem has been analyzed precisely in [23]. 

The interaction between DE and DM could be a major issue to be confronted in studying the physics of DE. However, 
due to the nature of these two components remain unknown, it will not be possible to derive the precise form of the 
interaction from first principles. One has to assume a specific coupling from the outset [24, 25] or determine it from 
phenomenological requirements [2G]. In view of the continuity equations, the interaction between DE and DM must 
be a function of the energy densities multiplied by a quantity with units of inverse of time. For the latter the obvious 
choice is the Hubble factor H. Thus, the interaction between DE and DM could be expressed phenomenologically in 
forms such as Q = Q{HpDM)[^^h ■">*'], Q = QiHpDE)[2S], Q = Q{H{pde + Pdj\/))[0]-[13] or most generally in the 
form Q = Q{HpoE , Hpdm) which leads Q — (5i HpoE + HpjjM from the first term in the power law expansion. 
Considering the couplings are terms in the lagragian which mix both DE and DM, one may further presume that they 
could be parameterized by some product of the densities of DE and DM, for example Q ~ ^PdePdm [U]- Besides 
these phenomenological descriptions of the interaction between DE and DM, recently there is an attempt to describe 
the coupling from the thermodynamical consideration[27, 28]. 

It is of great interest to investigate effects of the interaction between DE and DM on the universe evolution, especially 
the influence of different forms of the interaction on cosmological parameters. This is the main motivation of the 
present paper. Using the new 182 Gold type la supernova samples (SNIa), the baryon acoustic oscillation measurement 
from the Sloan Digital Sky Survey (SDSS-BAO) and the shift parameter of the Cosmic Microwave Background given 
by the three- year Wilkinson Microwave Anisotropy Probe observations (CMB shift), we are going to examine the 
efficiency of different coupling forms in solving the cosmic coincidence problem. In our study we will not specify any 
special model of DE. Considering recent accurate data analysis showing that the time varying DE gives a better fit 
than a cosmological constant and in particular, DE EoS can cross —1 around z = 0.2 from above to below[29], we 
will employ two commonly used parameterizations in our work, namely cui = loq + loiz/ {\ + z) = + wi(l — a) and 
ujii = ujQ + ioiz / {\ -\- zY = ujQ-\-{a — a?)uji. For comparison we will also examine the effects of the interaction between 
DE and DM for the constant EoS of DE. 



II. PHENOMENOLOGICAL MODELS DESCRIBING THE INTERACTION BETWEEN DE AND DM 

In the FRW universe, the Fricdmann equation reads 

^MlH"^ = Pr + ph + PDM+ PDE- (1) 
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Using the critical energy density = 3MpHQ, we can rewrite the Friedmann equation in the form 

= = nlil + z)' + + + pde{z)IpI + Pdm{z)IpI. (2) 

The energy density of the radiation is given by, 

Pr = O'b^cmb (3) 

where a\y is the Stefan-Boltzmann constant, T^mh = 2.12QK is the CMB temperature at present. pde{z) represents 
the DE density and pdm is the energy density of the cold DM. 
When DE and DM interact with each other, we have 

Pdm + iHpDM Q , (4) 
Pde + 2>H[1 + wde)pde = -Q , (5) 

respectively, where Q denotes the interaction term. Notice that the overall energy density of dark sectors pdm + Pde 
is conserved. 

From the equations above, phenomenological forms of the interaction between DE and DM must be a function of 
the energy densities multiplied by a quantity with units of inverse of time. The obvious choice of the quantity with 
units of inverse of time is the Hubble factor H. The energy densities could be chosen as density of either DE, DM, 
DE plus DM, the linear combination of DE and DM or the product of DE and DM. 

We will compare these different phenomenological expressions of the interaction between DE and DM and investigate 
their effects on cosmological evolutions and efficiencies in solving the coincidence problem. 

A. The coupling is proportional to the energy density of DE 

If the interaction term is proportional to the energy density of DE,we choose Q = SHpjjE and rewrite the continuity 
equations 4 and 5 

Pdm + SHpDM - SHpoE = (6) 

PDE + '-iHpDE{l+io) + SHpDE = (7) 

Employing the first commonly used parameterization of the EoS w/, we obtain a closed form solution for (7) 

PDE = p2,i;e-3-ie3-i/(l + ^)(l + ^)3(l+-o+-i)+5, (g) 

Insert this solution into Eq(6), we get 

PDM = (1 + zfpUi - Spl^il + zf re3../(i+.)-3..(i ^ (9) 



Choosing the second parameterization of the EoS loh, we have 

PDE = pl^e>^ + '2'-^/i^+-?-^^^/i^+-){l + ^)3(l+-o)+5^ (10) 

and 

PDM = pUA^ + zf-Sp^^il + zfei-^ /•%|.,/(l+.)^-3.,/(l+.)(i^^)^+3.o-lrf^. (11) 
B. The coupling is proportional to the energy density of DM 

Expressing the interaction between DE and DM in proportional to the DM energy density Q = SHpDM,^^ can 
rewrite the continuity equations as 

Pdm + "iHpDM - SHpoM = 0, (12) 

PDE + iHpDE{l+Uj) + SHpDM = 0. (13) 
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Employing as the parameterization of the EoS, from Eq(12) we obtain. 

PDAI^pUli^ + ^f^'- (14) 

Insert this pdm into Eq(13), we arrive at 

PDF = (l + z)3(i+-<'+-i)e3-i/(i+^)<5pV re"'"'/^'+'Hl + -2)"'*"''+"'^"'"'^^^ (15) 

Jo 

+p23£;(1 + z)3(i+"o+'^i)e3"i/(l+^)-3'^i . (16) 
Using ojji as the EoS, we find the evolution equation for DM as 

Pdai^pImC^ + z)'-'. (17) 
Then from Eq(13) we have the evolution of the DE reads 

PDE = (l + z)3(i+'^'')e"TTi+i(TTiF5p03^^ f {1 + zy^'^^-^-^e^'^TT^dz (18) 



"'0 

+ (l + z)3(i+"°'e"^+*uTiFe^"VDB- (19) 

C. The coupling is proportional to the total energy density of DE and DM 

We can also write the interaction between DE and DM in proportion to the total energy densities of DE and DM, 
so that the continuity equations read 

Pdm + iHpDM - SH{pdm + Pde) = 0, (20) 
PDE + 3Hpde{1 +oj) + SH{pDM + pde) = 0. (21) 

The above equations can be converted to 

- T^— (1 + ^)P£'£; - -4— (PiJB +Pi3m) = 0. (23) 
az I + z I + z 

D. The coupling is proportional to the linear combination of energy densities of DE and DM 

The phenomenological coupling between DE and DM could be expressed most generally in the form Q = 

Q{H Pde T H Pdm) which leads Q ~ SiHpDE + S2HPDM from the first term in the power law expansion. Thus 
the continuity equations read 

Pdm + 'iHpDM - H{5ipde + S2PDM) = 0, (24) 

PDE + iHpDE{l + W) + H{6ipDE + S2PDM) = 0, (25) 

which can be changed into 

dpDM _ 3 ^^^^ ^ 1 ^g^p^^ _|_ s^p^j^^-j = 0^ (26) 
az 1 + z I + z 

^ - t4-(1 + ""^PDE - T^iSlPDE + hPDM) = 0. (27) 

az \ + z \ + z 

For the last two cases of describing the interaction between DE and DM, it is not easy to write out analytic solutions of 
Pde and pdm when we substitute the parameterizations of EoS. We will completely count on the numerical calculation 
to investigate the time evolution of the DE and DM in these cases. 
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E. The coupling is proportional to the product of energy densities of DE and DM 



The coupling can be considered in terms of the lagrangian which mix both DE and DM, we can adopt the form 
of the interaction as the product of the densities of DE and DM [14], Q = \pdePdm- where A = kM~^ and M is 
considered as a mass parameter and k is a dimensionless constant. Taking Q.^ = Px/ Pco and p^o = •^HqAI^ / {Sir) = 

■ 0.81 • lQ~^^GeV'^, the dimensionless continuity equations for the DM and DE read 

dQpM 3 ^ ^ ^ 1 ^IpAiilpE _ Q 

dz 1 + 2 (1 + 2^) V^DAI + ^DE 

drtpE 3 , , 1 Q.DM^DE „ 

— 1 -r—0- + ^)ylDE-0— — -—==== = (28) 

dz l + z (1 + z) ^VIdm + ^DE 

where i5 = ^J^J{ST^\p].(^ Mp = ^2,/ (87r)fc(-^)(^)^/^, which is a dimensionless coupling constant now. 

So far we have listed all possible phenomenological descriptions of the interaction between DE and DM. In the next 
section, we will examine the viability of these descriptions by comparing with observations. 



III. OBSERVATIONAL CONSTRAINTS 



We wiU use the three-year WMAP (WMAP3) data[30], the SN la data [31] and the Baryon Acoustic Oscillation 
(BAO) measurement from the Sloan Digital Sky Survey [32] to study the property of the interaction between DE and 
DM. In [33], the authors showed that combining the shift parameters R and the angular scale la of the sound horizon 
at recombination appears to be a good approximation of the full WMAP3 data. The model independent constraints 
on R and la was later given [.iJ] by using the WMAP3 data, where the covariance matrix of the parameters i?, la and 
f2f,/i^ was also provided. In the following we will use the shift parameter R, the angular scale la of the sound horizon 
at recombination and their covariance matrix given in \'A] and employ the Monte-Carlo Markov Chain (MCMC) 
method to explore the parameter space. Our MCMC code is based on the publicly available package CMBeasy[38]. 

For the SNIa data, we will calculate = ^^"'"^^'^^ ^'^'^^ , where the extinction-corrected distance modulus 

p{z) — blogioldLiz) /Mpc] + 25, ai is the total uncertainty in the SN la data, and the luminosity distance is dL{z) = 
(1 + z) / Ho J^' dz'/E{z') and the nuisance parameter Hq is marginalized over with flat prior. 
For the SDSS data, we add the BAO parameter [30, 32], 

2/3 

= 0.469(0.95/0.98)-°-3^ ± 0.017, 
to x^- 

For WMAP3 data, we first add the shift parameter[34] 

, /"^'^ dz' 

R=V^l ^ = 1.70±0.03, 

to at z = 1089 ± 1. It was argued in [33, 34] that the combination of the shift parameter and the angular scale 
of the sound horizon at recombination gives much better constraints on cosmological parameters. Thus we add the 
angular scale of the sound horizon at recombination 

7ri?/^/no" 
7 — ' V — qn? 2-1-19 

QdzcjEiz)-''^-^^'-^' 

where the sound speed Cs = 1/ ■\/3(l + Rta), Rb — 315000ilh/i^(rcmfc/2.7if a is the scale factor, and rif,/i^ = 

n 099+0 00082 

U.UZZ_Q QQQg2. 

We use the MCMC method to explore the parameter space. We treat Hq as also a fitting parameter, and we 
impose a prior of Hq = 72 ± 8fcm/s/Mpc[35]. By using two parameterizations for the EoS of DE uji,ujii, we obtain 
the parameter space for different couplings respectively and we will show our results in the following sections. For 
comparisons we will also present the result for constant DE EoS. 



r-'^ dz' 



£;(0.35)i/3 0.35 7o E{z') 
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Figure 1: The likelihood for coupling parameters when the interaction is chosen as SHpDE 
A. The coupling is proportional to the energy density of DE 

Choosing the first parameterization of the DE EoS loi and expressing the coupling between DE and DM in propor- 
tional to the energy density of DE, the results are shown in Fig la and Fig 2a. The quantity of interest for analyzing 
the coincidence problem is the ratio r = pdm/ Pde and its evolution with time. The behavior of dr/dh\a deserves 
attention with regard to the coincidence problem. Clearly we see that the positive coupling obtained from the best fit 
data leads slower change of r which makes the coincidence problem less acute when compared with the case without 
interaction shown by the red dotted line in Fig 2a. When we choose DE EoS to be w// and constant, results are shown 
in Fig 2b, c. From the behaviors of the ratio r = pdm / PDE,ioT the constant DE EoS or cj//,this phenomenological 
interaction cannot help to alleviate the coincidence problem. Besides, the negative coupling shown in Fig lb for DE 
EoS ujjj leads negative Pde'^^ the early time. When DE EoS is constant pde is nearly zero even in the inflationary 
era, which makes it hard to count on this puE^o explain the infation. What's worse, the poM will be negative in the 
future. All these unphysical consequences are brought by the negative coupling from the fitting. 
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B. The coupling is proportional to the energy density of DM 

Using the first parameterization of the DE EoS and choosing the coupling between DE and DM in proportional to 
the energy density of DM, we show the likelihood in Fig 3a and the behaviors of cosmological parameters in Fig 4a. 
Comparing lines of dr/dlna for the interacting and nonintcracting cases, it seems that near the present time the 
interaction between DE and DM makes the change of r a bit slower than the nonintcracting case, which makes the 
coincidence problem less serious. However since the coupling is negative from the best fit. we encounter the negative 
Pde in the early epoch of the universe. This sounds unphysical[.>( i], although in some modified gravity theory negative 
value of PDE was allowed[37].The results for choosing DE EoS to be oju and constant are shown in Fig 4b and Fig 4c. It 
is clear from the behavior of r that the coincidence problem cannot be alleviated by choosing this phenomenological 
interaction form when DE EoS are of w// or constant. Moreover, negative poE very early epoch appears again. 
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Figure 2: These figures show clearly the behaviors of different dark sectors during the late time evolution of universe with 
best-fitted data for coupling SHpoE- The red dotted line denotes the noninteracting case. The left one corresponds to EoS 
[^/,the middle one is EoS un and the right one is constant EoS 
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Figure 3: The likelihood for coupling parameters when the interaction is chosen as SHpoM- 
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C. The coupling is proportional to the total energy density of DE and DM 

The results of the evolution of cosmological parameters for employing w/and expressing the interaction between 
DE and DM in proportional to the total energy density of DE and DM are shown in Fig 5a and Fig 6a. In Fig 6a 
we see that,rchanges a bit slower than that of the noninteracting case at very recent epoch, and this might help to 
alleviate the coincidence problem. However in the early epoch the negative coupling obtained from fitting makes pde 
negative again. Choosing DE EoS to be w//, besides the negative poE in the early time, poM will appear negative in 
the future. For constant DE EoS, p^E can aslo be negative in tlic early epoch. Furthermore from Fig 5b, c, we see 
that the coincidence problem has not been alleviated. 
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Figure 4: These figures show clearly the behaviors of different dark sectors during the late time evolution of universe with 
best-fitted data for the coupling SHpoM- The red dotted line denotes the noninteracting data.The left one corresponds to EoS 
tj/,the middle one is EoS ujii and the right one is constant EoS 




Figure 5: The likelihood for coupling parameters when the interaction is chosen as 5H{pdm + Pde)- 
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D. The coupling is proportional to the linear combination of energy densities of DE and DM 

If the coupling is in linear combination of energy densities of DE and DM ,we can see easily from Fig 9, the coincidence 
problem gets even worse , since r changes even faster than that of the noninteracting case.Furthermore,we find the 
negative pdm in the future for all the cases, which is unphysical. 
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Figure 6: These figures show clearly the behaviors of different dark sectors during the late time evolution of universe with 
best-fitted data for the coupling 5H{pdm + Pde)- The red dotted line denotes the noninteracting case. The left one corresponds 
to EoS Loi, the middle one is EoS uii and the right one is constant EoS 




Figure 7: The likelihood for coupling parameters when the interaction is H{SipDE + 52Pdm) 
.The left one shows the likelihood for parameter 5\ and the right one show the likelihood for parameter 82- The EoS is chosen 

as uji for these figures. 
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E. The coupling is proportional to the product of energy densities of DE and DM 

In Fig 11, we have shown the results by choosing the interaction as the product of energy densities of DE and 
DM w/. Choosing DE EoS uji,\t is clear that with this kind of interaction between DE and DM, there is a slower 
change of r as compared to the noninteracting case. This means that the period when energy densities of DE and 
DM are comparable is longer compared to the noninteracting case. Thus it is not so strange that we now live in 
the coincidence state of the universe. In this sense the coincidence problem is less acute when compared with the 
case without interaction. Adopting this kind of interaction form, we have not met unphysical problems as mentioned 
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Figure 8: The likelihood for coupling parameters when the interaction is H{5\pde + 52Pdm) 
.The left one shows the likelihood for parameter 5\ and the right one show the likelihood for parameter 82- The EoS is chosen 

as ujii for these figures. 
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Figure 9: These figures show clearly the behaviors of different dark sectors during the late time evolution of universe with best- 
fitted data for the coupling H{5ipDE + 52Pdm)- The red dotted line denotes the noninteracting case. The left one corresponds 
to EoS ci;/,the middle one is EoS uon and the right one is constant EoS 

above. Choosing DE EoS as wj/jWe see that in most period, r changes slower than that of the noninteracting case. 
Only very recently r changes a bit faster for the interacting case. This will not change the qualitative behavior of 
the period when energy densities of DE and DM are comparable, which is still longer compared to the noninteracting 
case. Thus the coincidence problem can become less serious. Choosing the interaction in the form of the product 
of densities of DE and DM, we have positive coupling from the fitting, which can help to avoid unphysical problems 
mentioned above. For constant DE EoS, we see that r changes a bit faster than that of the noninteracting case near 
the present epoch, which cannot help to alleviate the coincidence problem. However the positive coupling from fitting 
again avoids unphysical problems for the evolution of cosmological parameters. 
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Figure 10: The likelihood for coupling parameters when the interaction is chosen as XpoMpDE 
.The left one shows the EoS ui and the right one is uu 
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Figure 11: These figures show clearly the behaviors of different dark sectors during the late time evolution of universe with 
best-fitted data for the interaction XpoMpDE- The red dotted line denotes the noninteracting case. The left one corresponds to 
EoS oi/jthe middle one is EoS uon and the right one is constant EoS 



IV. CONCLUSIONS AND DISCUSSIONS 



In this paper we have examined the effects of all possible phenomenological interactions between DE and DM on 
cosmological parameters and their efficiency on solving the coincidence problem. We have worked with two simple 
parameterizations of the DE EoS and the constant DE EoS. Using observational data coming from the new 182 
Gold type la supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year 
Wilkinson Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan 
Digital Sky Survey, we have performed a statistical joint analysis of different forms of phenomenological interaction 
between DE and DM. We have found that for the time-dependent DE EoS the model of the interaction in product 
of densities of DE and DM is the best model to alleviate the coincidence problem. For other interaction forms, 
the consequences of DE and DM interaction on cosmological parameters are very sensitive to the DE EoS. For the 
DE EoS in the form of w/, choosing the DE and DM interaction in proportional to the energy density of DE, we 
found the solution to solve the coincidence problem. We observed that other forms of the interaction either bring 
unphysical problems or cannot help to alleviate the coincidence problem. For the DE EoS w//, we found that except 
the interaction in terms of the product of densities of DE and DM, other phenomenological descriptions of interaction 
between DE and DM fail to alleviate the coincidence problem. For the constant DE EoS, none of the phenomenological 
models of interaction can help to alleviate the coincidence problem. It would be interesting to extend our analysis to 
more general parameterizations for both the interaction and the DE EoS available in the literature. Furthermore it 
is also of interest to include more observational data to constrain the interaction between DE and DM and study its 
consequences on the cosmological parameters. 
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